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This  first  quarterly  report  summarizes  progress  under  the 
VELA  Network  and  Automatic  Processing  Research  Prgram,  Contract 
Numbei  F08606-75-C-0029,  during  the  period  1 December  1974  to  31 
March  1975.  Work  in  the  following  areas  is  summarized: 
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SECTION  I 

INTRODUCTION  AND  SUMMARY 


This  quarterly  report  summarizes  the  progress  made  dueing 
the  period  from  1 December  1974  to  31  March  1975  on  the  VELA  Network* 
Evaluation  and  Automatic  Processing  Research  program  being  carried  out  by 
Texas  Instruments  Incorporated  at  the  Seismic  Data  Analysis  Center  (SDAC) 
in  Alexandria,  Virginia.  -The  five  tasks  of  the  program  are: 

• Evaluation  of  the  Korean  Seismic  Research  Station  (KSRS) 

• Investigation  of  signal  detection  techniques 

• Investigation  of  signal  estimation  techniques 

• Investigation  of  discrimination  techniques 

• System  Engineering  Studies. 

Software  to  edit  and  reformat  KSRS  data  is  being  developed. 
When  complete,  it  will  allow  the  use  of  Texas  Instruments  long-  and  short- 
period  array  evaluation  packages  on  KSRS  short-period  data.  These  packages 
will  be  usee  to  evaluate  the  noise  and  signal  characteristics  at  the  Korean 
Array  during  the  remainder  of  the  contract. 

The  signal  detection  algorithms  previously  developed  have  been 
modified  to  establish  a fixed  alarm  rate.  A non-causal  quality  check  which 
eliminates  bad  channels  has  been  developed,  and  execution  time  has  been  re- 
duced. Future  plans  are  to  process  a large  group  of  data  to  determine  the 
best  detector  configuration, 

A study  of  the  adaptive  beamforming  technique  as  applied  to 
interfering  short-period  signals  has  been  completed,  and  a report  prepared. 
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Work  will  continue  on  a long -period  interfering  event  study,  and  on  the  develop- 
ment of  a distributed  signal  model. 

A data  base  and  software  for  the  short-period  regional  and  tele- 
seismic  discrimination  study  are  under  preparation.  Preliminary  evaluation 
of  some  first-zone  events  has  shown  that  P to  S phase  energy  ratio  is  an 
effective  discriminant.  The  Herrin-Tucker  network  magnitude  hias  model  has 
been  expanded  to  allow  for  non-uniformity  in  station  detection  thresholds,  and 
a report  has  been  written.  It  is  planned  to  spend  the  remainder  of  the  contract 
period  on  the  spectral  analysis  of  the  regional  and  teleseismic  data,  and  on  the 
time  and  frequency  domain  analysis  of  the  first-zone  data. 

Group  velocity  analysis  of  Rayleigh  and  Love  waveB  is  now 
possible  using  the  Interactive  Long  Period  Processing  SyBtem.  Evaluation 
of  the  performance  of  the  interactive  system  haB  begun,  and  work  is  complete 
on  interfacing  the  system  with  data  from  the  Seismic  Research  Observatories. 
In  the  future,  the  matched  filtering  package  will  be  modified  to  accomodate 
short-period  data. 
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SECTION  II 

PRELIMINARY  EVALUATION  OF  SHORT-PERIOD  KSRS  ARRAY 


A.  CURRENT  STATUS 

The  preliminary  evaluation  of  the  KSRS  short-period  array  will 
take  place  in  two  parts.  The  first  is  the  development  of  software  to  handle 
KSRS  data,  and  the  second  is  the  actual  analysis  of  the  data.  During  this  quar- 
ter progress  was  made  on  the  first  of  these  parts. 

An  edit  program  to  produce  KSRS  edit  tapes  compatible  with 
Texas  Instruments  short-period  NORSAR  software  has  been  developed  and  is 
being  currently  checked  out.  Required  modifications  to  the  NORSAR  software 
are  also  under  way.  NORSAR  short-period  programs  can  thus  be  used  to  an- 
alyze KSRS  data. 

Work  has  begun  on  a reformating  program  for  this  data.  When 
complete,  this  program  will  reformat  short-period  data  into  TI's  standard  long- 
period  format.  These  data  will  then  appear  on  tape  just  as  long-period  data  do. 
These  tapes  will  be  compatible  with  TI's  long-period  array  evaluation  package, 
after  some  modification  to  that  package.  Then  frequency-wavenumber  analysis, 
available  only  in  the  long-period  package, can  be  used  without  further  develop- 
ment. 

The  KSRS  data  available  at  th»  SDAC  covers  two  time  periods  - 
spring  of  1973  and  November,  1974.  Representative  data  samples  from  each 
period  have  been  plotted,  and  it  was  found  that  the  1973  data  are  probably  of 
unacceptable  low  quality.  The  data  from  1974  appear  to  be  normal. 
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B. 


FUTURE  PLANS 


In  the  immediate  future,  we  plan  to  complete  the  reformating 
program  and  the  modifications  to  the  existing  evaluation  package.  The  remain- 
der of  the  next  quarter  will  be  devoted  to  an  evaluation  of  the  noise  and  signal 
characteristics  at  KSRS,  and  the  preparation  of  a report. 


SECTION  III 

SIGNAL  DETECTION  TECHNIQUES 


A.  CURRENT  STATUS 

The  combined  Fisher  and  conventional  power  detector  algorithm 
developed  under  Contract  VT/4705  for  short-period  signals,  has  been  modified 
to  establish  a variable  detection  threshold,  which  maintains  the  alarm  rate  at 
a fixed,  predetermined  level.  Such  a threshold  was  found  to  be  essential  to  an 
optimum  detector  because  of  the  rapidly  varying  noise  levels  at  the  KSRS.  The 
quality  control  algorithm  has  been  made  non-causal,  so  that  it  drops  bad  chan- 
nels from  consideration  before  they  enter  the  beamforming  process.  A pre- 
liminary improvement  of  the  software  has  resulted  in  a reduction  of  the  ex- 
ecution time  to  about  two  thirds  of  the  real  duration  of  the  data  processed. 

B.  FUTURE  PLANS 

During  the  next  quarter  the  detectors  under  study  will  be  used 
to  process  seismic  events  from  the  Kuriles -Kamchatka  region,  if  an  adequate 
detection  log  for  this  region  can  be  obtained.  Otherwise,  data  from  the  whole 
of  Eurasia  will  be  used.  The  detector  configuration  with  the  best  probability 
of  detection  will  be  determined,  and  a report  written. 


Ill- 1 


SECTION  IV 
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SIGNAL  ESTIMATION  TASK 
A.  CURRENT  STATUS 

A study  of  the  adaptive  beamforming  technique  as  applied  to  in- 
terfering short-period  events  recorded  at  KSRS  has  been  completed  during  this 
quarter.  It  was  concluded  that  the  ABF  technique  produces  only  6 dB  gain  re- 
lative to  simple  beamforming  in  this  application,  because  of  poor  site-to-site 
signal  similarity  and  of  rapidly  changing  relative  channel  power.  This  study 
was  done  on  spring  1973  data,  and  it  might  change  significantly  if  high  quality 
data  were  used.  A report  has  been  submitted. 
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B.  FUTURE  PLANS 

During  the  next  quarter,  the  signals  estimation  task  will  concen- 
trate in  two  basic  areas.  The  first  is  a continuation  of  the  long-period  adaptive 
beamforming  interfering-event  study,  which  will  concentrate  on  the  effect  of 
changes  in  processing  parameters  such  as  azimuthal  and  time  separation  be- 
tween events.  The  second  area  is  Wiener  distributed  signal  model  adaptive 
filtering.  Effort  will  be  devoted  next  quarter  to  a theoretical  examination  of 
this  model  and  the  results  to  be  expected  from  the  application  of  the  adaptive 
beamforming  technique  to  such  signals. 
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SECTION  V 
DISCRIMINATION 


A.  CURRENT  STATUS 

1.  Part  of  the  discrimination  task  includes  the  investigation  of  the 
power  of  short-period  NORSAR  source  spectrum  measurements  as  discrimi- 
nants. 

A data  base  for  this  study  is  in  preparation,  which  includes 
earthquakes  and  explosions  at  regional  and  teleseismic  distances.  A program 
to  filter  traces  and  compute  spectra  by  the  maximum  entropy  method  is  near- 
ing completion. 

2.  The  purpose  of  the  first-zone  discrimination  study  is  to  deter- 
mine the  capability  of  a Eurasian  seismic  station  to  discriminate  among  events 
at  epicentral  distances  of  less  than  2000  kilometers  using  phase  energy  ratios. 

A frequency-time -power  analysis  was  performed  on  near-field  events  record- 
ed at  NORSAR.  In  every  case  the  ratio  of  the  S phase  energy  to  P phase  energy 
discriminated  correctly  between  earthquakes  and  presumed  explosions.  The 
data  base  for  this  study  is  being  expanded. 


3.  A method  to  extend  the  Herrin-Tucker  network  magnitude  bias 

model  to  allow  for  non-uniformity  in  the  station  detection  thresholds  has  been 
developed.  This  method  applies  maximum  likelihood  estimation  to  determine 
event  network  magnitude,  given  the  individual  magnitudes  at  the  detecting  sta- 
tions and  the  measured  noise  levels  at  the  non-detecting  stations.  It  was  found 
that  an  upward  bias  of  from  a half  to  one  magnitude  unit  (mb)  exists  in  the  or- 
dinary magnitude  estimate  for  events  near  the  threshold  of  a 10  element  net- 
work. The  maximum  likelihood  method  removes  this  bias.  A report  has  been 
submitted  on  this  subject.  This  completes  the  contractural  work  on  this  task. 
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B. 


FUTURE  PLANS 


1.  NORSAR  short-period  source  spectrum  measurements  in  the 
next  quarter  will  include  corner  frequency,  roll  off,  and  long-period  displace- 
ment. Noise  and  system  response  corrections  are  major  problems  in  comput- 
ing accurate  spectra,  while  absorbtion  and  diffraction  may  create  problems  in 
their  interpretation.  These  problems  will  be  dealt  with  in  the  coming  quarter. 
Significant  deviations  of  the  corrected  explosion  sources  from  accepted  scale 
laws  for  earthquakes  will  be  taken  as  discriminants. 

2.  During  the  coming  quarter  we  will  complete  the  editing  of  our 
additional  first-zone  long-period  NORSAR  data  and  do  a visual  analysis  of  both 
short-  and  long -period  data  to  measure  travel  times,  amplitudes,  and  magni- 
tudes. Moving  window  spectral  analysis  will  be  performed  on  all  data  and  the 
energy  ratios  between  various  phases  will  be  estimated.  The  effectiveness  of 
these  quantities  as  discriminants  will  be  measured  and  reported  on. 
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SECTION  VI 

SYSTEM  ENGINEERING  STUDIES 
A.  CURRENT  STA'  ,S 

The  Inte  active  Long  Period  Processing  System  (ILPPS)  for 
the  PDP-15  computer  has  been  developed  further  toward  becoming  an  opera- 
tionally flexible  interactive  matched  filtering  package.  An  option  to  perform 
group  velocity  analysis  of  Rayleigh  and  JLove  waves  has  been  integrated  with 
the  previously  developed  system,  and  an  option  to  generate  a hard  copy  of  the 
information  displayed  on  the  CRT  screen  has  been  implemented.  It  is  possible 
with  this  option  to  generate  a hard  copy  using  either  the  CALCOMP  plotter  or 
the  VARIAN  electrostatic  printer /plotter. 

Evaluation  of  the  performance  of  the  interactive  system  has 
been  iniated  to  determine  its  operational  practicality  and  what  changes,  if 
any,  are  required.  Both  the  spectral  analysis  package  for  group  velocity  de- 
termination and  the  interactive  complex  cepstrum  technique  for  extraction  of 
later  phases  and  multiple  events  are  currently  being  evaluated. 


Work  is  complete  on  interfacing  the  interactive  system  with 
data  tapes  generated  by  the  Seismic  Research  Observatories  (SRO). 

The  capabilities  now  exist  to: 

• Input  a specified  short-  or  long-period  event  or  noise  sample 
directly  from  an  SRO  field  data  tape  to  a disk  file. 

• For  long-period  data,  the  full  signrl  analysis  package  develop- 
ed in  the  ILPPS  system  is  available  to  SRO  data,  once  the  above 
disk  files  have  been  generated. 
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For  short-period  dafa,  the  capability  exists  to  perform  time 
window  selection,  bandpass  filtering,  and  parameter  measure- 
ments along  the  same  lines  as  in  the  ILPPS  system. 


B.  FUTURE  PLANS 


Evaluation  of  the  interactive  techniques  will  continue,  with  em- 
phases on  the  expected  performance  in  an  operational  surveillance  system. 
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